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Abstract 

Spin current, i.e. the flow of spin angular momentum or magnetic moment, has recently 
attracted much attention as the promising alternative for charge current with better en¬ 
ergy efficiency. Genuine spin current is generally carried by the spin wave (propagating 
spin precession) in insulating ferromagnets, and should hold the chiral symmetry when 
it propagates along the spin direction. Here, we experimentally demonstrate that such a 
spin wave spin current (SWSC) shows nonreciprocal propagation characters in a chiral- 
lattice ferromagnet. This phenomenon originates from the interference of chirality between 
the SWSC and crystal-lattice, which is mediated by the relativistic spin-orbit interaction. 
The present finding enables the design of perfect spin current diode, and highlights the 
importance of the chiral aspect in SWSC. 
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Electron is a particle characterized by the charge and spin degree of freedom. 
In contrast to the charge current accompanied by Joule heat loss, spin current 
can be dissipationless and potentially minimize the energy consumption associated 
with the information processing'll - . Spin current is generally carried by the spin- 

polarized conduction electrons in metallic system, as well as the spin wave in insulat¬ 
ing system js, 6|. In particular, the latter spin wave spin current (SWSC) possesses 
many advantages over the former one, since it can avoid the simultaneous flow of 
charge current and has much longer propagation length. From the viewpoint of the 
symmetry, there are two types of SWSC for ferromagnets, depending on the direc¬ 
tional relationship between the carried magnetic moment Mo (parallel to the external 
magnetic held H) and the wave vector k. In case of M 0 || k (Fig. HJ4), the SWSC 
doesn’t have any mirror plane or space-inversion center, and belongs to the chiral 
symmetry 7], 8]. In contrast, the SWSC with M 0 J_ k configuration (Fig. [Tp) has 
the polar symmetry with the polar axis normal to both M 0 and k. For each case, 
the reversal of k gives the SWSC with opposite chirality or polarity (Fig. CD B and 
D). 

The above analysis predicts that the SWSCs propagating along the positive and 
negative direction can show different propagation characters, when placed in the chi¬ 
ral or polar environment depending on the symmetry of SWSC. For example, the 
surface (or interface) is always characterized by the structural polarization normal 
to the surface. When both Mq(\\ H) and k are confined within the surface plane 
keeping the M 0 J_ k relationship, the polar axis of SWSC (Fig. CD C and D) becomes 
parallel or antiparallel to the polarization of the surface and thus the asymmet¬ 
ric spin wave propagation between ±k can be expected. Such a surface-induced 
spin wave nonreciprocity has first been predicted by Damon and Eshbach more 
than 50 years ago[9], and then verified by various experimental techniques such as 
spin wave spectroscopy jlO, H], spin-polarized electron energy loss spectroscopy 12], 
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thermography [l3|, and Brillouin light scattering[l4]. In particular, the k —> 0 mode 
mediated by the magnetic dipole-dipole interaction in this configuration is called 
magnetostatic surface wave (representing that it can propagate only at the surface 
of the sample), and is known for its unidirectional propagation character | q, QQ. 
Nevertheless, such a surface-driven nonreciprocity mostly cancels out as the entire 
sample, since the sign of polarity is generally opposite between the top and bottom 
surfaces. 

In contrast, the interplay between the SWSC and the chiral medium has hardly 
been investigated. One reason is that the coexistence of crystallographic chirality 
and ferromagnetism is very rare in the real materials. Recently, however, some 
metallic ferromagnets with the chiral B20 crystallographic lattice (such as MnSi, 
Fei-zCo^Si, and FeGe) have been found to host magnetic skyrmions, i.e. nanometer¬ 
sized vortex-like swirling spin texture with particle nature lfF19|. Skyrmions are 
now attracting much attention as the potential information carrier for the high- 
density magnetic storage device 19,l20|, and the above hireling has promoted further 
search of similar chiral-lattice ferromagnetic materials. Soon after their discovery in 
the B20 materials, it has been reported that a chiral-lattice ferromagnetic insulator 
Cu 2 0Se0 3 (with cubic space group P2i3) can also host magnetic skyrmions, for 


the narrow temperature region just below T c ~ 59 K 


21-24], Notably, this material 


offers an ideal opportunity to investigate the property of SWSC under the chiral 
environment. In this work, we have examined the nonreciprocal nature of spin wave 
propagation for this chiral-lattice magnet in terms of spin wave spectroscopy. 

The basic concept of the spin wave spectroscopy[25|, [26] as well as the employed 
device structure is summarized in Fig. [U E and F. A pair of Au coplanar waveguides 
(ports 1 and 2) were fabricated on the oxidized silicon substrate, and the plate¬ 
shaped single crystal of Cu20Se03 was placed across them. Here, the chirality of 
each Cu20Se03 crystal (D or L) is checked by measuring the natural optical activity 
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at light wavelength 1310 nm in advance. When the oscillating electric current I v of 
gigahertz frequency v is injected into one of the waveguides, I v generates oscillating 
magnetic field H v and excites spin wave (i.e. coherent magnetization oscillation M v ) 
in the Cu 20 Se 03 sample. The propagating spin wave causes an additional magnetic 
flux on the waveguides, and induces the oscillating electric voltage V v following the 
Faraday’s law. By measuring the spectrum of complex inductance L nm {y) as defined 
by V” = (with m and n representing the port numbers used for the 

excitation and detection, respectively) with the vector network analyzer (VNA), we 
can evaluate both magnitude and phase of propagating spin wave. The spin wave 
contribution to the inductance spectrum A L nm [y) = L nm {y) — U^iy) is derived 
by the subtraction of the common background U^Jy) from the raw data L nm (u). 
Here, L nm (y) taken at H — 2650 Oe is adopted as L^(u), where the magnetic 
resonance is absent within our target frequency range from 2GHz to 7GHz. The 
wave number k of excited spin wave is determined by the spatial periodicity A(= 
12yum) of the waveguide pattern and the associated current density I u (x) [25|, 26]. 
Its Fourier transform \I u {k)\ 2 has the main peak at k p = 0.50/im -1 with the full 
width at half maximum (FWHM) of 5k = 0.37/mi -1 as plotted in Fig. 0P, satisfying 
the relationship k p ~ 2tt/X. To investigate the property of SWSC of the chiral 
symmetry (Fig. [U A and B), the H || k || [001] configuration is always adopted 
here. In this setup, the k —>■ 0 mode is called magnetostatic backward volume wave, 
which propagates through the entire volume of the sample with the negative group 


velocity^, [ill, Il5j. In the centrosymmetric materials, this mode should not show any 


nonreciprocal propagation nature. 

First, we have investigated the nature of SWSC in the uniform collinear ferro¬ 
magnetic state with saturated magnetization. Figure [2l A and B indicate the real 
and imaginary part of ALn and AL 2 i spectra measured at +740 Oe, i.e. in the 
collinear ferromagnetic state, for the D-chirality of the Cu 2 0 Se 03 crystal. The self- 
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inductance A Ln represents the efficiency of the local spin wave excitation, and the 
ferromagnetic resonance characterized by Lorentzian shape of spectrum can be iden¬ 
tified at 3.2 GHz. In contrast, the mutual-inductance AL 21 reflects the propagation 
character of spin wave between the two waveguides, and the finite oscillating sig¬ 
nal can be detected around the same resonance frequency. Hereafter, we focus on 
the comparison between AL 21 and AL 12 , each of which stands for the propagating 
spin wave characterized by the wave vector +k and —k, respectively. To interpret 
the data more intuitively, the spectra of |AL nm | and 0 as defined with A L nm = 
Re[AL nm ] ± i Im[AL nm ] = | AL nm | exp[i0] are plotted in Fig. [2J C and D. |AL 2 i| 
and |ALi 2 | express the magnitude of spin wave after the propagation along the pos¬ 
itive and negative directions, and both spectra show a peak structure. Notably, the 
peak frequency u p as well as the peak intensity |AL^ m | are clearly different between 
±k. On the other hand, 0 represents the phase delay of spin wave for a transmission 
between the two waveguides separated by the distance d(— 20/iin), and will satisfy 
the relationship 0 = kd when a single spin wave mode is assumed|25l. 26]. This means 
that the 0 spectrum directly reflects the spin wave dispersion relationship, and its 
slope gives the group velocity v g = ^ = 27rd(^) _1 . The clear difference in the 0 
slope between the ones derived from AL 2 i and AL 12 indicates that v g of spin wave 
is not equal between ±k. The above results establish that the propagation character 
of spin wave in this configuration is nonreciprocal, from both aspects of magnitude 
and group velocity. 

To reveal the origin of nonreciprocity, the measurements of Im[AZ/ 2 i] and Im[ALi 2 ] 
were performed with various combinations of the magnetic field direction (H = ±740 
Oe) and crystallographic chirality (D and L) for Cu 2 0Se0 3 (Fig. [2] E to H). Fig¬ 
ure [2]E shows the data for the D-crystal at ±740 Oe. The spectra of Im[AL 2 i] and 
Im[A±i 2 ] are characterized by the signal oscillating with different period and magni¬ 
tude, in agreement with the feature observed for | AL nm | and 0. For the opposite sign 
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of applied H (Fig. |2F), the spectral shapes for Im[AL 2 i] and Im[ALi 2 ] (i.e. the sign 
of nonreciprocity) are reversed. Likewise, the employment of opposite chirality of 
crystal (i.e. L-crystal) also reverses the sign of nonreciprocity (Fig. [2] G and H). In 
Fig. [2l I to L, the symmetrically expected sign of nonreciprocity for each experimen¬ 
tal configuration is summarized. In general, the SWSC is expressed as the product 
of magnetic moment M 0 (|| H) and wave vector k (Fig. [TJ A and B). This means 
that //-reversal has the same effect as fc-reversal, and thus results in the opposite 
sign of nonreciprocity. On the other hand, the application of space-inversion oper¬ 
ation to the system reverses both the crystallographic chirality and the ^-direction 
keeping H and Mo unchanged, indicating that the sign of nonreciprocity should be 
opposite between L-crystal and D-crystal. The above symmetry-based analysis is in 
agreement with the experimental results, which proves that the observed spin wave 
nonreciprocity originates from the chiral nature of crystal lattice. 

Next, we investigated the magnetic field dependence of nonreciprocity. In Fig. [3l 
A and B, the peak frequency v p for |AL 2 i| and |ALi 2 | (defined as u p (+k) and u p (—k)), 
as well as the difference between them (A u p = u p (+k) — u p (—k)), are plotted as a 


function of H. Cu 2 0 Se 03 is known to host the helical spin order for H = 0 


21 


23 |, 


while it is replaced with the uniform collinear ferromagnetic order for H > 600 Oe. 
The magnetic resonance frequency is gradually suppressed by H in the helical spin 


state, and then shows //-linear increase in the collinear 


behaviors are consistent with the previous reports 24], 


27 


erromagnetic state. These 


28]. Notably, the magni¬ 


tude of nonreciprocity is essentially dependent on the underlying magnetic structure. 
While the relatively large shift of resonance frequency Av p ~ 0.05 GHz between ±k 
is always observed for the collinear ferromagnetic state, such a nonreciprocity sud¬ 
denly vanishes upon the transition into the helical spin state. Similar behavior is 
also observed for the peak intensity |AL^ m | (Fig. [3]F) and group velocity v& (Fig. 
W) deduced at u p . For all these properties, clear nonreciprocity is observed only in 
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the collinear ferromagnetic state. Note that u pi and v p basically reflect the 

frequency, magnitude, and group velocity of spin wave for \k\ = k p . At maximum, 
the frequency shift Au p = 0.07 GHz and the change in v p (|AL£ m |) up to 25 % (40%) 
can be obtained between ±k. Their sign of nonreciprocity is confirmed to reverse for 
//-reversal. 

To clarify the microscopic origin of observed nonreciprocity, we attempt 


tirnate the spin wave dispersion for this material. According to Ref. 29, 30], the 


o es- 


magnetic Hamiltonian for the ferromagnets with chiral cubic lattice symmetry under 
the continuum approximation can be written as 


n = 


VS) 2 -DS • [V x S] 


K 




dr ., 


( 1 ) 


where J, ZZ, and K describe the magnitude of ferromagnetic exchange, 
Dzyaloshinskii-Moriya (DM), and cubic anisotropy term, respectively. S is dimen¬ 
sionless parameter representing the vector spin density. 7 , /i 0 , h = 2irh, and V 0 
are gyromagnetic ratio, vacuum magnetic permeability, Planck constant, and the 
volume of formula unit cell of Cu 20 Se 03 , respectively. For the H || k || [001] config¬ 
uration, the spin wave dispersion u(k) for the uniform collinear ferromagnetic state 


is described as 


29| 


u = ^[2 DSk + JSk 2 + 2 KS 3 } + —fi 0 H. 
h 27r 


( 2 ) 


In the real sample, this dispersion is further modified by the additional contribution 
of the magnetic dipole-dipole interaction, especially for the k —>■ 0 region. When the 
infinitely wide plate-shaped sample with the thickness l is assumed and H || k || [001] 
lies along the in-plane direction, Eq. [2] can be rewritten asjll, 15, 30] 


2DSVnk 1 


v = 


f -\ (JSV 0 k 2 + 2KV 0 S 3 + ~fh/i 0 H)(JSV 0 k 2 + 2 KV 0 S 3 + -jh t i 0 H + 7 


h 


( 3 ) 


(1 — e-l fc F). 

[ 
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with M s being the saturation magnetization. In Fig. [H D and E, the spin wave 
dispersion calculated based on Eq. [3] with the material parameters estimated for 
CuoOSeOa is plotted. Equation [3] can be approximated by Eq. [2] except for the 
k —> 0 region, and gives parabolic dispersion with its minimum at k = —D/J. As 
k approaches zero, however, the contribution of magnetic dipole-dipole interaction 
gradually increases the spin wave frequency. It causes the negative group velocity 
for the k —> 0 region, and this mode can be considered as a kind of magnetostatic 


backward volume wave 
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, 3 - 


Note that the first and second term in Eq. [3] are odd 
and even functions of k, respectively, and thus only the former one proportional to 
Dk can contribute to the spin wave nonreciprocity. This suggests that the observed 
nonreciprocity directly comes from the DM interaction, whose sign and magnitude 
reflect the chirality of the underlying crystallographic lattice through the relativistic 
spin-orbit interaction. 

Experimentally, the above spin wave dispersion relationship can be partly repro¬ 
duced by analyzing the </> spectrum (Fig. [2D). Given that the frequency v — v p 
corresponds to the wave number k = k p and the relationship (j) = kd holds 25, [26], 
the dispersion relationship can be determined by k = [4>(v) — (j){y p )\/d + k p . In Fig. 
0]A, the spin wave dispersions v{k) for the collinear ferromagnetic state deduced from 
the AL 21 and AL 12 spectra in Fig. [2j C and D are plotted, each of which corresponds 
to the one for positive and negative k, respectively. For both cases, the concave-up 
dispersions with negative slopes are obtained, consistent with the prediction of Eq. 
[3] for the k —> 0 region (Fig. l4E)fll. 3 . 30], The dispersion curves for positive 
and negative k show considerable deviation from each other, and their frequency 
shift Av(\k\) = v(+\k\) — v(—\k\) is plotted in Fig. [4]B. Av is found to be almost 
proportional to \k\, which is consistent with the relationship Av = 4DSVo\k\/h ex¬ 
pected from Eq. [3] The observed slope of Av gives D ~ 5.5 x 10~ 4 J/m 2 , which 
roughly agrees with D ~ 3.4 x 10~ 4 J/m 2 estimated from //-dependence of mag- 












netic resonance frequency 27]. These results firmly confirm that the observed spin 
wave nonreciprocity stems from the DM interaction associated with the chiral nature 
of crystallographic lattice. Note that for the helical spin state under the magnetic 
Hamiltonian given by Eq. [U it has been proposed that the Brillouin zone is folded 
back with the helical spin modulation period due to the expansion of magnetic unit 


cell 


29 


31]. Such a folding back of magnon branch should extinguish the asymmetry 


between ±k, which explains the observed disappearance of spin wave nonreciprocity 
in the helical spin state. 

The above relationship A v oc D\k\ obtained for the ferromagnetic state sug¬ 
gests that Av linearly increases for larger \k\. Since the linewidth of the resonance 
peak in A L nm spectrum is mainly determined by the ^-distribution of waveguide 
for the present situation (see Supporting Online Material), the employment^ 
waveguide pattern with shorter wavelength and/or repeated meander shape 


o' 

the 

25, 

26] 


(i.e. enhancement of k p and suppression of 5k) will completely resolve the frequency 
overlap of spin wave signal between ±k. This means that for the given resonance 
frequency the SWSC can propagate along only one direction, and not along the 
opposite direction at all. 

Recently, the relevance of DM interaction has also been discussed for the case of 


32 


34], through the spin-polarized 


14] experiments. Their reported 


surface/interface-driven spin wave nonreciprocity 
electron energy loss 12] and Brillouin light scattering 
magnitudes of frequency shift between ±k are Av/\k\ = 50 MHz//imfor Pt/Co/Ni 
him [nj] and 240 MHZ//Z 1 I 1” 1 for Fe double layer on W(110) 12], while the presently 
observed 130 MHz//im -1 for bulk Cu 2 0Se0 3 is comparable with these systems. Note 
that the nonreciprocity in C^OSeOa is for the volume spin wave and originates from 
the chiral lattice symmetry of the bulk crystal itself, unlike the conventional case 
of the surface/interface-driven nonreciprocity for the surface spin wave. Since the 
nonreciprocal volume spin wave can avoid the cancelation of nonreciprocity between 
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the top and bottom surface inherent to the latter situation, our present finding 
offers a simple and promising route for the realization of perfect spin current diode 
with the 100% efficiency of rectification. It may also find unique spin caloritronic 


applications js, 


135 


, such as magnetically-tunable unidirectional heat conveyer 


Q. 


From a broader perspective, any (quasi-)particle flow along the magnetic held 
direction should have the chiral symmetry{3], and will show the similar nonreciprocal 
propagation character in the chiral-lattice compound. This_has been experimentally 
confirmed for the light 


36 


37] and conduction electron 


38 


, while their reported 


magnitude of nonreciprocity is generally very small. Our present observation of clear 
spin wave (or magnon) nonreciprocity in the chiral-lattice compound highlights the 
importance of the chiral aspect in SWSC. 
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SUPPORTING ONLINE MATERIAL 


m 


Experimental Details 

Single crystals of Cu20Se03 are grown by chemical vapor transport method 


4lj . They are cut into the plate-like shape with widest faces parallel to the (110) 


plane, and polished with diamond slurry and colloidal silica. D- and L-chirality of 
crystals are distinguished by measuring the sign of natural optical activity at light 
wavelength 1310 nm. Every piece of crystal shows optical rotation angle ±16°/mm, 
and its single-domain nature is confirmed by the observation under polarized-light 
microscope. 


Our procedure of the measurement and data analysis associa 
wave spectroscopy basically follows the method proposed in Ref. 


ea 


23 , 


with the spin 


26]. A pair of 


coplanar waveguide patterns (ports 1 and 2) consisting of Au 195 nm / Ti 5 nm are 
deposited on a thermally oxidized silicon substrate using photo lithography and elec¬ 
tron beam evaporation technique. The rectangular shape of CuoOSeOa crystals with 
typical size from 60 /.ini x 20 /jrn x 2 /jrn to 60 /jrn x 10 /jrn x 1 /jrn are extracted 
from the original bulk crystal pieces by focused ion beam micro-sampling technique, 
and placed across the waveguides with W deposition at both edges of the crystal. 
This device is put into the probe station equipped with GM refrigerator and hori¬ 
zontal electromagnet, and connected with a vector network analyzer (VNA) through 
the coaxial cable and GSG (ground-signal-ground) microprobe. The calibration is 
performed using Short-Open-Load-Throurgh coplanar standards. The spectrum of 
^-parameter ( S nm (u , H) with m and n representing the port numbers used for the 
excitation and detection, respectively) is measured by VNA at various magnitudes 
of external magnetic field, and converted into the impedance spectrum Z nm (is,H) 
assuming the characteristic impedance Zq = 5011 42], The one at H re f = 2650 Oe is 
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considered as the background, where no magnetic resonance appears within the tar¬ 
get frequency range (2 GHz < v < 7 GHz). From the subtraction of two impedance 
spectra, the spin wave contribution to the inductance spectrum AL nrn (v, H) is de¬ 
rived as A L nm (u, H ) = [Z nm (u, H ) — Z nm (y , H ve f)\/(i27ru). The input power into the 
waveguide is —15 dBm within the linear response regime, which is confirmed by mea¬ 
suring the power dependence of spectrum around the magnetic resonance frequency. 
In this work, the H || k || [001] configuration is always adopted. 


DESIGN OF COPLANER WAVEGUIDE AND LINEWIDTH OF MAG¬ 
NETIC RESONANCE 


In Fig. [5]A, the coplanar waveguide pattern employed in this study is illustrated. 
The wavelength A and propagation gap d are 12 /mi and 20 /mi, respectively. Each 
waveguide consists of one signal line at the center and two ground lines at the both 
sides, which is terminated with a short circuit. When it is connected to the VNA 
through the GSG microprobe, the input current density for the signal and ground 
line is Jq and -Jq/ 2, respectively. By taking the Fourier transform for the spatial 


distribution o 
estimated 


25 


current density I u (x), the wavenumber distribution \I u (k)\ 2 can be 
26] as shown in Fig. [5j3- The main peak at k p = 0.50pm -1 satisfies 
the relationship k p ~ 27T/A, and its full width at half maximum (FWHM) is 5k = 
0.37 pm -1 . The higher order peaks are also found for the larger k region, and the 
second largest one is at k = 1.47pm -1 with the amplitude 7 times smaller than 
that for the main peak. To simplify the discussion, we analyzed our A L nm spectra 
assuming that the contribution from the main peak centered at k p is dominant. 

Such a wavenumber distribution in waveguides (Fig. [5j3) directly affects the 
linewidth of ferromagnetic resonance. In Im[ALn] spectrum, the FWHM 5v for the 
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( 4 ) 


resonance peak at frequency u p can be given as[26] 


v p a ■ 5k 

5u = -2- -1- 2v p a, 

Zn 

with a representing the intrinsic Gilbert damping parameter. In case of the present 
Cu 20 Se 03 specimen at 740 Oe, 5v = 0.42 GHz is obtained from the Im[ALn] 
spectrum of Fig. 2A in the main text. Considering the corresponding averaged spin 
wave group velocity v v g = 6.1 km/s taken from Fig. 3D in the main text, the first 
term in Eq. [4] gives ~ 0.36 GHz. This means that 5v mostly reflects the wave 
number distribution associated with the waveguide pattern. By using v p = 3.2 GHz, 
we obtain the relatively small damping parameter a ~ 0.01, which is consistent with 
the previous report [27]. This allows us to estimate the decay length of propagating 
spin wave Id = v p /(2-k au p ) = 30/irn 26]. 


Material properties 


Our tar get material, Cu 20 Se 03 has the chiral cubic crystal lattice with space 
group P2i3 43, 44J. It contains two distinctive magnetic Cu 2+ (S = 1/2) sites with 
the ratio of 3 : 1, and three-up one-down type of local ferrimagnetic spin arrange¬ 


ment has been reported below magnetic ordering temperature T c 


59 K 


44 , m ]• 


This material hosts helical spin order under zero magnetic field, where spins rotate 


within a plane normal to the magnetic modulation vector q 21-23]. Application 


of magnetic field aligns q parallel to H and turns the spin texture into the conical 
one. Above the critical magnitude of magnetic held H Cl magnetization is saturated 
and collinear ferrimagnetic state is stabilized. In the present study, the character 
of spin wave in this collinear ferrimagnetic state is mainly investigated. To simplify 
the analysis, we adopted the continuum approximation and treated this compound 
as the ferromagnetic system. 
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The material parameters included in Eq. 3 in the main text can be estimated 
so as to reproduce the //-dependence of magnetic resonance frequency v p in the 
ferromagnetic state (Fig. [6]) . In this process, several additional confinements are 
imposed 29]. The helical spin modulation period A h (~ 62 nm) 22, .23] and the 
corresponding magnetic wave number Q = 2ir/\ h in the ground state is given as 


and the critical magnetic field 

7 h 


Q = -D/J , 

; (~ 0.063 T) satishes 

D2S 2 KS>. 


( 5 ) 


J 


( 6 ) 


The saturation magnetization is M s = h'yS/V o = 0.46/tb/Cu 2+ at 30 K, with Vo ~ 89 
A 3 being the volume of formula unit cell of Cu 20 Se 03 . From these restrictions, we 
obtain D = 3.4 x 10 -4 J/m 2 , J = 3.4 x 10 -12 J/m, K = 6.9 x 10 3 J/m 3 , 7/27r = 29 
GHz T _1 , HqM s = 0.12 T, and S = 0.44. These values are used to calculate the spin 
wave dispersion in Fig. 4, D and E in the main text and //-dependence of v p in Fig. 
[6l Similar values have also been reported in Ref. 0. 

Note that Cu 20 Se 03 also hosts the skyrmion spin state for narrow temperature 
region from 56 K to 59 K, just below T c \ 21, 2l|, |23]. Since the magnitude of prop¬ 
agating spin wave rapidly decays as temperature approaches T c , the experimental 
investigation of spin wave character in the skyrmion state remains the future chal¬ 
lenge. 


Symmetry of Spin Current 

Spin current is generally characterized by the combination of magnetic moment 
Mq and wave vector k. Here, Mq is an axial vector, and odd for time-reversal and 
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even for space-inversion. In contrast, A: is a polar vector, and odd for both time- 
reversal and space-inversion. Tims, the spin current expressed as the product of k 
and Mo is even for time-reversal and odd for space-inversion. 

Figure [71 A and B summarize the compatible symmetry elements for the spin 
current with the k || Md and k _L Mo configurations. In case of k || Mo, the spin 
current sustains a rotation axis along the k direction and 2' (two-fold rotation followed 
by time-reversal) axis normal to it (Fig. [7jA) . Since no mirror plane or space-inversion 
center is present, this belongs to the chiral (and not polar) symmetry. In contrast, 
the spin current with the k _L Mo configuration has a mirror plane (m) normal to Md, 
m' (mirror reflection followed by time-reversal) plane normal to k , and 2'-axis normal 
to both k and Mq (Fig. [7]B). This belongs to the polar (and not chiral) symmetry 
with the polar axis normal to both k and Mq. 

For each configuration, the application of space-inversion operation reverses k 
(but not Mo), as well as the associated sign of the chirality or polarity in the spin 
current. 

Figures 
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FIG. 1: (A)-(D) Spin waves for the uniform collinear ferromagnetic state, characterized 
by various combinations of the wave vector k and the uniform magnetization component 
Mo, which can be considered as the flow of magnetic moment (i.e. spin current). Such 
a spin wave spin current belongs to chiral (polar) symmetry for the k || Mo (k _L Mo) 
configuration, and the reversal of k gives opposite sign of chirality (polarity). The red and 
gray arrows represent the directions of local magnetization and its precession, respectively. 
(E) The optical microscope image and (F) schematic illustration of the device structure 
used for spin wave spectroscopy. In (E), the directions for positive sign of H and k are also 
indicated. 


18 























Cu 2 0Se0 3 /-/|| A: || [001] 30 K 
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Frequency (GFIz) 



FIG. 2: (A)-(H) Spin wave contribution to inductance spectrum A L nm =Re[AL nm ] + i 
Im[AL nm ] = \AL nm \ exp[i0], measured for the D- or L-chirality of Cr^OSeOs single crys¬ 
tal with the FI || /c || [001] configuration at 30 K. All the data are taken at the uniform 
collinear ferromagnetic state. (A) and (B) Real and imaginary part of self inductance 
ALn and mutual inductance AL 21 measured for the D-crystal at H = +740 Oe. For the 
same configuration, (C) magnitude |AL nm | and (D) phase </> of AL 21 and AL 12 are also 
plotted. (E)-(H) Imaginary part of AL 21 and AL 12 , measured with various combinations 
of magnetic field direction (H = ±740 Oe) and crystallographic chirality (D or L). Note 
that the deviation of the overall signal magnitude and resonance frequency between D- 
and L-crystal is due to the slight difference in their sample size and associated demagne¬ 
tizing field. The corresponding experimental configurations as well as the expected sign of 
nonreciprocity are summarized in (I)-(L). Here, the spin wave characterized by the wave 
vector +k (—k) contributes to AL 21 (AL 12 ), and the solid and dashed arrows represent 
the different propagation characters. 
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Magnetic Field (Oe) 


FIG. 3: Magnetic field dependence of spin wave nonreciprocity between AL 21 and AL 12 
(i.e. +k and —k), measured for the D-crystal of Cu 20 Se 03 with the H \\ k \\ [001] 
configuration at 30 K. (A) and (B) indicate the magnetic resonance frequencies fp giving 
the peak value of |AL nm |, and their difference between ±k (i.e. Au p = v p (+k) — z/ p (— k)), 
respectively. In (C) and (D), the corresponding peak value \ALnm\ and the group velocity 
Vg at the frequency u p are also plotted. 
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D-crystal H||/c||[001] 
30 K +740 Oe 



FIG. 4: (A) Spin wave dispersion for the D-crystal of Cu 20 Se 03 at H = +740 Oe (i.e. 
collinear ferromagnetic state), experimentally deduced by analyzing the 4> spectrum in Fig. 
ED. The one obtained from AL 21 (AL 12 ) corresponds to positive (negative) A;- value, and the 
frequency difference between ±k (i.e. Azz(|fc|) = v(+\k\)—i'(—\k\)) is also plotted in (B). (C) 
Wave number distribution of excitation current I u (k), obtained by the Fourier transform of 
the wave guide pattern. (D) and (E) Spin wave dispersions for the collinear ferromagnetic 
state calculated based on Eq. [3] or Eq. El the latter of which ignores the effect of magnetic 
dipole-dipole interaction. The black dashed line represents the contribution of the second 
term in Eq. [3j Here, the assumed material parameters are l = 2/xm, D = 3.4 x IIP 4 J/m 2 , 
J = 3.4 x 10 -12 J/m, K = 6.9 x 10 3 J/ 111 3 , 'y/2ir = 29 GHz T _1 , /J>oM s = 0.12 T, = 
0.074 T, E 0 = 89 A 3 , and S = 0.44. 
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FIG. 5: (A) Schematic illustration of a pair of coplanar waveguides used for the spin wave 
spectroscopy. Each waveguide consists of one signal (S) line and two ground (G) lines. The 
associated current density distribution as well as length scale (in the unit of /im) are also 
shown. (B) The calculated wavenumber distribution of excitation current. 
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D-crystal H||/c||[001] 30 K 



Magnetic Field (Oe) 


FIG. 6: Magnetic field dependence of resonance frequency v p in the AL 21 spectrum, mea¬ 
sured for the D-crystal of Cu 20 Se 03 at 30 K. The experimental data is taken from Fig. 
3A in the main text, and the theoretical fit is by Eq. 3 in the main text with k = k p . 
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FIG. 7: Spin current characterized by the combination of magnetic moment Mq and wave 
vector k , with (A) k || Mq and (B) k _L Mq configurations. The compatible symmetry 
elements are also indicated. 


23 









